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Abstract 

Topological quantum phases of matter are characterized by an intimate relationship between the 
Hamiltonian dynamics away from the edges and the appearance of bound states localized at the 
edges of the system. Elucidating this correspondence in the continuum formulation of topological 
phases, even in the simplest case of a one-dimensional system, touches upon fundamental concepts 
and methods in quantum mechanics that are not commonly discussed in textbooks, in particular 
the self-adjoint extensions of a Hermitian operator. We show how such topological bound states 
can be derived in a prototypical one-dimensional system. Along the way, we provide a pedagogical 
exposition of the self-adjoint extension method as well as the role of symmetries in correctly for¬ 
mulating the continuum, held-theory description of topological matter with boundaries. Moreover, 
we show that self-adjoint extensions can be characterized generally in terms of a conserved local 
current associated with the self-adjoint operator. 
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I. INTRODUCTION 


Phases of matter are operationally distinguished by their relevant physical observables. In 
quantum mechanics, these observables correspond to linear, self-adjoint, and gauge-invariant 
operators acting on a properly specified Hilbert space, whose elements represent the physical 
states of the system.^ In this formalism, measurable values of an observable are identihed 
with the eigenvalues of the corresponding operator, and since the result of a measurement 
is associated with a real number it is sensible to impose the condition of self-adjointness to 
guarantee real-valued eigenvalues. When the Hilbert state space is inhnite-dimensional, some 
operators may be unbounded and, therefore, need to be treated with special care; in particu¬ 
lar, determining the proper set of boundary conditions that establishes their self-adjointness 
becomes a crucial step in their dehnition.^ Specihcally, for an unbounded operator to be 
self-adjoint, it is necessary but not sufficient for it to be Hermitian (symmetric). There¬ 
fore, a naive approach to analyzing the eigenvalues and eigenvectors of Hermitian operators 
without properly specifying the self-adjoint boundary conditions can lead to unexpected 
paradoxes,^ resulting in wrong physical conclusions. 

These considerations become quite important for the study of topological phases of mat¬ 
ter,® where the properties of the physical system near boundaries are closely linked with 
those of the bulk. In view of the increased interest in topological phases, it is a goal of 
this work to illustrate the correct implementation of the theory of self-adjoint extensions 
of operators to the continuum, or held-theory, description of topological quantum matter. 
Fortunately, there is already a well-developed theory of self-adjoint extensions of Hermitian 
operators® however, its application to, and consequences for, the continuum description 
of topological phases remain obscure. We believe this is in part due to the mathematically 
abstract nature of this theory. Thus, we also aim to provide a formulation of the theory of 
self-adjoint extensions that is more natural for physicists. 

We will illustrate this method in the context of a concrete and simple model of a topologi¬ 
cal phase. From a pedagogical standpoint, topological phases provide a very natural setting 
for learning and using self-adjoint extensions. It is a topic that brings together several 
strands of modern physics and mathematics: many-body quantum systems; symmetries of 
a physical system; topology; and functional analysis. We present an accessible exposition of 
these subjects in the context of a simple model of a ferromagnetic insulator with spin-orbit 
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interactions. We also clarify the role of symmetry and boundary conditions in self-adjoint 
extensions of an operator. In particular, we show that the notion of bulk-boundary cor¬ 
respondence, i.e., a relationship between the topological properties of the bulk phase and 
the character of its boundary states (see, e.g.. Chap. 6 of Ref. |3l), depends on the choice 
of self-adjoint extension. Finally, we provide a new and intuitive physical interpretation 
of self-adjoint extensions in terms of a conserved current associated with the self-adjoint 
operator. 

The paper is organized as follows. In Sec. [IT] we briefly explain the notions of a topolog¬ 
ical state of matter and its concomitant boundary states. To illustrate these concepts, in 


Sec. Ill we introduce a model of a one-dimensional lattice that exhibits normal and topo¬ 


logical phases over a range of parameters. In Sec. IV we present a pedagogical derivation 
of the continuum, held-theory description of this model with boundaries, and point out the 
difficulty of dealing with unbounded operators. In Sec. Mwe review von Neumann’s theory 


of self-adjoint extensions of a Hermitian operator. In Sec. VI we use this theory to obtain 


the self-adjoint extensions of the continuum model Hamiltonian of Sec. [IV] and its topolog¬ 
ical bound states. We also discuss the effect of a symmetry operation on the self-adjoint 


extension. In Sec. VII we recast the self-adjoint extensions of arbitrary Hermitian operators 


in terms of a conserved, spatially local current. We conclude in Sec. VHI Technical details 
are presented in four appendices. 


II. TOPOLOGICAL QUANTUM PHASES AND BOUNDARY STATES 

The study of topological phases of quantum systems has flourished in recent years, thanks 
to theoretical and experimental discoveries of several families of such phases in artihcial 
nanostructures as well as bulk crystals.® In contrast to the phases of matter characterized 
by spontaneously broken symmetries, such as ferromagnets vs. paramagnets or solids vs. 
liquids,^ topological phases of electronic matter are characterized by the nontrivial topology 
of their electrons’ wave functions. In particular, a topological phase cannot be distinguished 
from a normal phase by probing its bulk electronic properties locally. For example, a topo¬ 
logical insulator does not allow the passage of electric current through it just like a normal 
insulator; however, as one changes the momentum of the electrons, in a topological insulator, 
their wave functions exhibit a nontrivial twist. This twist cannot be undone by changing the 
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parameters of the system unless the system is brought to a critical phase at which it becomes 
a metalPThis characterization also implies that at the boundary between a topological phase 
and a normal phase, the electronic motion must undergo a drastic change. For instance, the 
electronic states at the boundary between a topological and a normal insulator are metallic. 
In other words, one expects that at the boundaries of the system in a topological phase there 
are bound states, protected by symmetry, that emerge because of the nontrivial topology of 
the bulk. This statement is commonly referred to as the bulk-boundary correspondence. 

Though crystals are described microscopically in terms of their lattice structure, it is 
often more convenient to describe these systems in a continuum, held-theory, representa¬ 
tion. Such a description is relevant for the long-distance behavior of the system, which is 
what determines its quantum phases. It is also a more general description since different 
lattice systems can end up having the same continuum description and, thus, the same 
long-distance physics. As is generically the case in quantum systems, continuum descrip¬ 
tions involve unbounded operators such as the linear momentum. Usually this does not 
present a major obstacle since one can regularize such operators by restricting their action 
to square-integrable wave functions that vanish sufficiently fast at inhnity. However, in or¬ 
der to illustrate the bulk-boundary correspondence in a topological phase, one must study 
the properties of the system near the boundaries. In this case, one needs to regularize the 
unbounded operators differently. Intuitively, this regularization can be seen to involve the 
boundary conditions of the electronic wave functions. Mathematically, it requires the notion 
of self-adjoint extensions of unbounded Hermitian operators 

In order to illustrate these concepts, in the next section we introduce a simple lattice 
model in one spatial dimension, which exhibits a normal and a topological insulating phase. 
We derive the continuum Hamiltonian description of this model and study the spectral 
properties of the corresponding held theory. We shall see that correctly identifying the 
phases of the system in this continuum formulation requires the proper identihcation of 
boundary conditions related to the self-adjointness of the Hamiltonian operator. 


III. A TOPOLOGICAL QUANTUM WIRE 

Consider a system of spin-1/2 fermions moving along a line with a periodic array of po¬ 
tential wells, for example electrons moving through a chain of N ions, subject to a magnetic 
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field and an internal spin-orbit interaction. When the separation between the ions a is large 
enongh so that only the qnantnm tnnneling between nearest ions is appreciable, one can 
model this system as a discrete chain with the Hamiltonian,® 
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( 1 ) 

Here, “-fh.c.” means one mnst add the Hermitian conjngate of all the previons terms. We 
are nsing a second-qnantized notation of the creation and annihilation operators in the Fock 
space, and whose action is to create and remove, respectively, an electron at ion 
site r G {1, 2, • • •, N} with spin s G {t, i} along the direction of the magnetic held. These 
operators obey the anti-commutation relations c^^CsV' + = 0 and c^As'r' + — 

Sss'Srr', where 6 is the Kronecker delta. The operator product cAsr counts how many 
electrons (0 or 1) of spin s are at site r. The operator product displaces an electron 

of spin s from site r to one of spin s' at site r -|-1. In this way, we can see that the parameter 
fi in the model rehects the Zeeman energy splitting in the magnetic held; w is an energy 
scale related to the strength of the spin-dependent quantum tunneling between two nearest 
ions; and A is related to the strength of the spin-orbit interaction. (Note the absence of 
a spin-independent tunneling term, and the fact that the sum in Eq. ([^ does not include 
terms with Cj,jy+i and When the number of electrons is N, the system may represent 

a ferromagnetic insulator with spin-orbit coupling. For simplicity, we shall take all the 
parameters to be real and w and A to be positive numbers. 

Note that this model has the following property: upon hipping the spin c-fr i—)■ c^r and 
A hnds H I—>• —H. Denoting the operation with the unitary map S = (5^)“^, 

we can express this as SHA = —H. This means that if we hip the spins in any energy 
eigenstate \E) with energy E, H\E) = E\E), we obtain an eigenstate with energy —E, since 


H{S\E)) = -{ShAs\E) = -SH\E) = -E{S\E)). (2) 


Therefore the spectrum of H is symmetric around E = 0. 

This quantum wire is a prototypical system that displays diherent quantum phases upon 
changes in the parameters of the Hamiltonian that are distinguished by their topological 
properties. This can be seen more easily if we rotate the spin basis to a direction orthogonal 
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FIG. 1. The single-particle energy spectrum of the quantum wire shows the quantum phases of 
the system. Here, N = 150 and X/w = 5.5. In the topological insulator phase, \fi/2w\ < 1, a pair 
of zero-energy bound states is found localized at the two edges of the lattice. They spilt away in 
the normal insulator phase. 


to the magnetic field. In the new basis, the annihilation operators are 

£ hj-j. -(- 

O^r ~ 


C^r 


V2 • 

After some straightforward algebra, we obtain the Hamiltonian in the new basis as 

N r 




r=l 


+ ( ^ + ^) K+lhr + ^ 


-b h.c. 


( 3 ) 

( 4 ) 

( 5 ) 


Now we can easily see that when the parameters are at the “sweet spot” /i = 0 and w = A/2, 
the Hamiltonian reduces to H = X Therefore, the two operators and 

(and their Hermitian conjugates) do not appear in the Hamiltonian and [H, 6'|-i] = 
[H, b^isf] = 0. This means that the bound states created by and at the two ends of 
the wire are degenerate eigenstates of the Hamiltonian (with zero energy). 

One can show that even away from the sweet spot, there are two bound eigenstates of the 
Hamiltonian as long as |/x| < 2w. The energy of these bound states is exponentially small 
in the system size, Ei, ~ Here a is the lattice spacing between nearest-neighbor ions 

and ^ is a length scale that depends on system parameters. In the thermodynamic limit, 
N —)■ oo, we recover two asymptotically exact zero-energy bound states. The existence or 
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absence of these bound states marks the topological or trivial phases of the system. The 
existence of these phases is illustrated in Fig. which shows the single-particle energy 
spectrum of the lattice Hamiltonian as /i is varied. 

Instead of studying the phases on the lattice, we focus on the thermodynamic limit in 
the following and derive a continuum Hamiltonian. This approach has the advantage of 
simplifying the solutions since it is often easier to solve continuous differential equations 
rather than discrete difference equations. Also, many different lattice models can have the 
same continuum description in the thermodynamic limit. Thus, by keeping only the relevant 
long-distance, universal features of the system in the continuum limit, we will gain a more 
general perspective on the potential phases of the system. 

IV. THE CONTINUUM LIMIT 


We shall take the continuum limit of the lattice model as the lattice constant a —)■ 0 
while the length L = Na is hxed. We leave open the choice of whether the length is hnite, 
semi-inhnite, or inhnite. In this limit, the discrete site index r is turned into a continuous 
position X. A simple way to deduce the relation between fermion operators Cgr in the lattice 
and the held operator iii fhe continuum is to demand the proper anti-commutation 

relations 


Ts(x)Ts/(a;') T5/(a:')Ts(a:) = 0, (6a) 

= 6ss'S{x - x'), (6b) 


where 6(x — x') is the Dirac delta function.^ Thus, the fermion held Ts(a;) must have the 
dimension of (length)“^A_ Dehning x = ar, x' = ar', and \im.a^o{5rr'/a) = ^(a: — a;'), we see 
that the proper dehnition of continuum operators has the form fl/s(a;) = \im.a^Q{csr/\/a). In 
the limit a —?• 0 we can expand 


-^Csr+i '^six + a) = Ts(a:) a^Ts(x) -h O(a^). (7) 

^/a ax 

Also in this limit —)■ J dx. Replacing these limits in the Hamiltonian, Eq. ([^, we hnd 

several terms that, to hrst order in a, vanish upon integration. In particular, terms of the 
form 


d 


d u 


aw 


Tj(a;)—T^(a;) dx-1-h.c. = aw / — T^(x)T^(x) 


dx 


dx 


dx 


( 8 ) 
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vanish in the limit a —)■ 0 for hnite w. After some straightforward algebra, we hnd the 
continuum Hamiltonian 

H ^ H^= ! ii\x)Hii{x)dx, (9) 


where 


'h(a;) = 




n = 


m V:r 

ax 


—v-^ —m 

ax 


( 10 ) 


with m = 2w + fi and v = aX. Note that in this limit we keep v hnite. Using the Pauli 
matrices 


— 


0 1\ fO -i 

u = 

10/ Vi 0 

the matrix Ti could be written in a compact form as 

H = tvay -—h mcr^. 
ax 


(Tz = 


1 0 
0 -1 


(11) 


( 12 ) 


We note that this is the Hamiltonian leading to the one-dimensional Dirac equation.^l The 


spin-hip operation of the lattice Hamiltonian, S of Sec. HI, is now implemented by the Pauli 
matrix S = ax and ST-iS"^ = —"H. 

The matrix 1-L is an operator in the position basis. Its eigenvalues E and eigenstate^ 
4>e{x) = (j)iE{x)y are found by solving the linear system of diherential equations 


'H(j)E{x) = E(()e{x). 


Using these solutions, we can dehne the energy creation operators 


(13) 




(j)^E{x)'^\{x) + (j)iE{x)'^l{x) 


dx. 


(14) 


As we show in Appendix the are “ladder operators” that satisfy the commutation 
relations [Hc^ In this way, the full spectrum of the continuum Hamiltonian can 

be obtained by solving the hrst-quantized eigenvalue equation ( [l^ . We will concentrate on 
studying the properties of di in the rest of this paper. 


The careful reader may have noticed that the energy E in Eq. (13) can be unboundedly 
large. For instance, consider a solution such as 


(f)E{x) = e 


ikx I 


E—m 

ivk 


(15) 







with eigenvalue E = E\/v‘^k‘^ + The energy E can be an indefinitely large number for 
large k. Thus Ti and the continuum Hamiltonian He, unlike the lattice Hamiltonian, are 
unbounded operators. More surprisingly, even though the Hamiltonian is formally Hermi- 
tian, in a finite geometry (a; < oo and/or x > —oo) the energy eigenvalues can be imaginary 
numbers. We can see this by replacing k with m in the above solution. Such a solution 
would be acceptable since e~'^^ need not diverge in a finite geometry for a proper choice of 
sgn(/i;). When \vk\ > \m\, the energy is a purely imaginary number. 

The presence of imaginary eigenvalues is problematic for the correspondence of the Hamil¬ 
tonian to observable energy. Formally, it signals the loss of self-adjointness of the Hamilto¬ 
nian operator. To restore self-adjointness, the Hamiltonian needs to be properly defined on 
a subset of the Hilbert space such that no imaginary eigenvalues exist in the spectrum. In 
the following two sections we will introduce the mathematical technique to deal with this 
problem for a general unbounded operator, and find the appropriate subset that resolves 
this issue for our continuum Hamiltonian. 

The following exercise shows that there is a second continuum field operator that one can 
define as $*( 0 :) = \mia^o{—^Ycsrly/a. 


Exercise 1. By expanding the field operators in the limit a —)■ 0, show that 
the continuum Hamiltonian for $ = ($'[■ is 




d 

+ (/i - 2w)az 


<F(a;) dx. 


(16) 


Thus, as the analysis in the following sections shows, the continuum limit in terms of <F 
becomes important for the topological phases of the system when ji — 2w changes sign. For 
simplicity, we shall assume 2w > fi in the following. 


V. THE SELF-ADJOINT EXTENSIONS OF A HERMITIAN OPERATOR 

The state of a quantum system is specified by a vector 0 in a Hilbert space The 
Hilbert space is equipped with an inner product (010) that is a comp lex-valued, positive- 
definite, sesquilinear function on ^ x ^ that maps any two states 0, 0 G ^ to a complex 
number (0|0) = (0|0)*. The inner product defines a norm for the Hilbert space, ||0|| = 
For example, given an interval [Li,L 2 ] C K, the Hilbert space ^ = ^‘^[Li,L 2 \ is 
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the set of all square-integrable functions with the inner-product ('?/’|0) = 'll)*{x)(f){x) dx, 

i.e., those functions with ||0|| < oo. 

Physical observables of a quantum system are specihed by linear self-adjoint operators.^ 
For an operator A to be self-adjoint, it primarily must be Hermitian (or symmetric). An 
operator A is dehned to be Hermitian if 

iijlAcj)) = {Am (17) 


for all vectors "0 and 0 in its domain. In many textbooks on quantum mechanics, no 
distinction is made between Hermitian and self-adjoint operators. Indeed, one need not 
make such a distinction for bounded operators. (A bounded operator B is one for which 
there exists a positive number M such that ||i?0|| < Af||0|| for all vectors 0 in the Hilbert 
space.) However, the distinction between Hermitian and self-adjoint operators becomes 
crucial for unbounded operators. This is so because for unbounded operators the designation 
of the domain of the operator A, S^{A), which is a subspace of the full Hilbert space 
&{A) C generally becomes an important part of its dehnition. 

For Hermitian operators, the connection between unboundedness and a restricted domain 
is quite general. In fact, according to the Hellinger-Toeplitz theorem,II^a Hermitian operator 
B with the domain ^{B) = is bounded. It is a direct corollary of this theorem that an 
unbounded and Hermitian operator cannot be dehned on the whole Hilbert space. Let us 
illustrate this point by an example. Consider the position operator, X, on the Hilbert space 


.if^[0,L] dehned as (X0)(a:) = a;0(a;) for 0 G .if^[0,L]. Equation (17) is satished since 


(0|X0) = / 'll)*{x)x(j){x) dx = / (a;0(a;))*0(a;) dx = (X0I0) 


(18) 


Jo Jo 

Thus, the position operator is a Hermitian operator. In this case, X is a bounded operator 
since for any 0 G .if^[0,L], ||X0|| < M||0|| with M = L. Now, if instead we consider the 
whole real line, R, the position operator dehned on the Hilbert space continues to be 

Hermitian, but is no longer bounded. Now, consider 0 G .if^([R) such that 0(x) = x/{x‘^ + 1). 
It is clear that (X0)(a;) = x^ j (x^-l-l) is not square-integrable and therefore X0 is ill-dehned. 
Therefore, in this case X cannot be dehned on the whole Hilbert space and one must hnd 
an appropriate domain for X, such as ^(X) = {0 G : ||X0|| < cx)}. This means 

0(x) —)■ 0 for X ^ oo faster than x~^^‘^. As we shall discuss in detail below, restrictions on 
the domain of an operator are commonly expressed in terms of boundary conditions on the 
state vectors. 
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In order to define a self-adjoint operator properly, we first need to define the adjoint of 
an operator and its domain. For an operator A : ^(A) —)■ and 0 G ^(A), the adjoint A'^ 
and its domain of states xjj G ^(A'^) are dehned by 




(19) 


(This dehnes nniquely only if the operator A is “densely dehned;” see Appendix for 
a discussion.) A self-adjoint operator A is then dehned as one that is Hermitian and has 
the same domain as its adjoint, ^(A^) = ^(A). 

As an example, consider the momentum operator, P, on the Hilbert space =2^^[0,L], 
dehned as (P0)(a;) = —id(p/dx (we use units such that Planck’s constant h= 1). Since the 
derivative of a square-integrable function can be arbitrarily large, P is unbounded. For its 
domain we take 3>{P) = {0 G .if^[0,L] : is absolutely continuous and(;/)(0) = 4’{L) = 

0}.!^ (One can show that with this dehnition, P is densely dehned; see Appendix]^) Simple 


algebra shows that Eq. (19) holds true if we have 

{P^'ip){x) = -i^i){x), 

r{LmL) = r{o)m- 


(20a) 

(20b) 


The hrst equation shows that the momentum operator is a Hermitian operator. However, 


the second equation shows that V’(O) and ip{L) could be nonzero and still Eq. (19) would 
be satished. That is, .^(Pf) = {■0 G =2f^[0,L] : 'ip{x) is absolutely continuous} D ^(P). 
Therefore, P is not self-adjoint. Some puzzling consequences of the lack of self-adjointness 
for P are discussed in Ref. El An important consequence for our purposes is that the domain 
of the adjoint operator Pi is so unconstrained that it admits imaginary eigenvalues. That 
is, one can hnd states V’± ^ ^(P^) such that PV± = for > 0. Indeed, one can 

easily show that the two states 


^P+{x) = e-^", ^P-{x) = (21) 

are such eigenstates with equal norm, ||t/’+|| = 

As seen in this example, the domain of the Hermitian operator A is no larger than the 
domain of its adjoint. Ah This means that in order to make the operator A self-adjoint 
one must extend its domain in a way that shrinks the domain of the adjoint to match it, 
and eliminate the imaginary eigenvalues from the spectrum of the adjoint operator. If the 
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procedure is successful, we find a self-adjoint extension of the Hermitian operator A. This 
procedure was made systematic by von Neumann by a method known as the dehciency 
indices method. 

In this method, given a Hermitian, densely dehned, and closed operator A : ^(A) —)• 

(see Appendix for the dehnitions), one hrst hnds the adjoint operator and its domain 
and then solves for the imaginary eigenstates of the adjoint operator. These solutions 
dehne the two deficiency subspaces = {-0^ G A^ip± = The dimensions 

of these subspaces, n± = dim(,j^), are the dehciency indices of A. 

Deficiency indices theorem (Von Neumann) For a Hermitian, densely 

dehned, and closed operator A, 

1. If n+ = n_ = 0, A is self-adjoint. 

2. If n+ 7 ^ n_, A is not self-adjoint and it has no self-adjoint extensions. 

3. If n+ = n_ = n > 0, A has inhnitely many self-adjoint extensions An 

parametrized by n x n unitary matrix maps U : —)■ For a choice 

of U, S>{A) c S>{Au) C ^(At), 

^(^w) = {'fi = +'Ip++ : (j) e ^(A), V’+ e je+}, (22) 

Aufi = AV = A(j) + iri'fi+ — if] Ufi+. (23) 


Note that even though rj appears explicitly in the dehciency subspaces, the family of self- 
adjoint extensions is independent of the parameter rj^ So, sometimes the choice rj = 1 is 
made. However, we will keep r] for dimensional purposes. 

Considering the example of the momentum operator P on .^^[0,L], we see that n_|_ = 
n_ = 1. Therefore there is a one-parameter family of self-adjoint extensions Pu = Pe with 


W : -0+ H-)■ and given in Eq. (|2l|), such that 


Pe^ = ( 24 ) 

^{Pe) = {^ = 0 -F c(0+ + e'V-) : 0(0) = 0(T) = 0, c G C}. (25) 


Note that G ^{Pe) satishes a new boundary condition, fi{L) = e*40(O), with = 
+ e^%l {I+ e-^^e^%. 
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While this method is systematic, it is rather abstract. In Sec. VII we will introduce a 
more intuitive way to hnd the proper domain of self-adjoint extensions in terms of boundary 
conditions and the local current associated with the corresponding self-adjoint operator. 


VI. SELF-ADJOINT EXTENSIONS AND TOPOLOGICAL BOUND STATES 


We will now apply von Neumann’s dehciency indices method to obtain self-adjoint ex¬ 


tensions of the continuum Hamiltonian "H, Eq. (12), in semi-inhnite and hnite-size wire 


geometries. In each case, we will show that the appropriately extended domains admit ad¬ 
ditional bound states in the topological phase, and we will relate their properties to the 
boundary conditions at the ends of the wire. We will also see that demanding that the spec¬ 


tral symmetry, given by the spin-flip operation S of Sec. |IV[ be preserved by the self-adjoint 
extension forces the bound states to have zero energy. 


A. Semi-infinite wire 

The simplest geometry with a boundary is the semi-inhnite wire, x G [0, oo), with just 
one physical boundary. In order to apply the dehciency theorem, we choose the domain 

^{n) = !</>= G cx)), 0^(0) ^ 

where s is the spin. (This domain is dense; the proof is similar to the example given in 
Appendix!^. With this choice, the adjoint operator has the domain : 'rps{x) G 

cxd)}. Now we can hnd all possible extensions T-iu such that ^{Hu) = 

Next, we solve for = ±iri'ip±. It is easy to see that 

V>±W = e-W-/”( Y (27) 

where |e|e*^ = m — irj. (Notice that 6 ^ 0,7r and i|Y’+|| = IIY’-IIO So the dehciency indices 
are = n_ = n = 1. The domain of the extended operator Tiu = 'He is parametrized by 
the one-parameter mapping U : i—>■ 

^{He) = -F c(Y’+ + e'V-) : Y>(0) = 0, c G C}. (28) 
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The domain is equivalently characterized by a new boundary condition. To see 

this, note that for an arbitrary state -0 G ^{He) we have 


^(0) 


^t(0)\ l + e‘® A 


(29) 


Thus 


A 


^ cos(g/2) 
cos(6'/2 — (5) 


(30a) 

(30b) 


Since we have extended the domain of "H, we may expect the spectrum of "Hg to depend 
on the choice of 9 or, equivalently, A. Indeed, one can easily check that for mA > 0 there is 
always a bound state. 


'ijjhi.x) = exp 


^ 2A X/a\ 
u A2 + l^yl^J’ 


(31) 


with energy eigenvalue 

^ A2- 1 

= aTTT’"' 

Note that —\m\ < Eb < \m\, i.e., it is within the energy gap of the bulk states. 

As the careful reader may infer from the dependence of the spectrum on the choice of self- 
adjoint extension, other physical properties may also depend on this choice. An important 
case is the effect of a symmetry operation. For example, what happens to the spectral 
symmetry given by the spin-flip operation S = ax of Sec. [ryl ? It is easy to see that under 
the operation S the boundary condition characterizing S>{'H 0 ) is mapped as A i—>■ A' = 1/A. 
Therefore, S is not in general a symmetry of the spectrum and the spectrum may not be 
symmetric around zero. Indeed, the existence of a single bound state with Eb 0 is possible 
only if the spectrum is not symmetric. The spectral symmetry is restored for A' = A = ±1; 
for A = sgn(m) the spectrum has a zero-energy bound state, Eb = 0. Therefore, whether or 
not a symmetry of the bulk lattice is a symmetry of the continuum description depends on 
the choice of self-adjoint extension. 
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B. Finite wire 


We now consider a more realistic wire geometry: a wire with a finite length L. We choose 


the domain of the Hamiltonian in Eq. (12) to be 


»(W) = e L], .#..(0) = ML) = 0}, 


(33) 


where s is the spin. (One can show that this is a dense domain similar to previous cases.) 
With this choice, the domain of the adjoint operator is 


»(«*) = W:V>.e^"10,L|}. 


(34) 


The deficiency indices are determined by solving the imaginary eigenvalue equation 'H^'ijj± = 
±ir]il>±. We now find two independent solutions for each sign, 





i’±\x) 



(35) 


with lele**^ = m — if] as before. (Note also that ||'0+^|| = ||'0-^ll = IIV'+^II = 11'*/’-^11-) The 
deficiency spaces = {V’± = c^±4>±* + G C} are two dimensional; the 

deficiency indices are n+ = = n = 2. Consequently, the extended domain is parametrized 

by 2 X 2 unitary matrix maps 


U ; 




(36) 


The domain of the self-adjoint extension is 


= {iJ = 4) +'Ip++ WV’+ : </>s(0) = 4)s(,L) = 0, tlL+e X+}. (37) 

What boundary conditions characterize To find the answer, we note that for 

any y,'0 G we must have 


(xl'Hw'0) = i'Huxl'ip)- (38) 

In particular, this must be true if we choose y = -0+^and y = separately. 

By partial integration, it is easy to see that 




+iv \4j^±\xf ay4j{x) 


a = 1,2. 


(39) 
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The boundary terms are straightforward to calculate and yield 




exp 


Ca-1 

2 



yi(L)+C„e’’'V.t(L)] 


^ ±d“’. 


( 40 ) 


with C,a = From Eq. (38), collecting the boundary terms, we hnd the matrix equation 

where ^+["0] = and ^-[ip] = Since we know that the 

family of self-adjoint extensions is the same for all choices of t], we can rewrite this for a 
particular choice of t]; namely we shall take the limiting case 77 —)■ oo (i.e., |e| —)■ 00 and 


e z) to hnd 


u 


+#t(0)^ _ / V’i(0)+#t(0) 

V + #t(^) / 


(41) 


The boundary condition can also be written in a form that appears independent of lA. 
This can be done by noting that the boundary conditions relate the two vectors and 
by a unitary matrix u. Equivalently, for two states G ^(Tiu) the inner product 
between the two vectors [- 0 ] and [ 0 ] dehned separately for -0 and 0 must be preserved, 
i.e., ^+['0]l.^+[0] = After some algebra this yields 


i’iWML) - i’l(L)ML) = y(o)^i(o) - ■q(0).#.t(o). 


(42) 


This condition is quadratic in the wavefunctions and so is not immediately useful in solving 
the differential equations needed to hnd the spectrum. However, it is equivalent to the 


linear boundary conditions in Eq. (41). We note that it reduces to the boundary condition 


we found for the semi-inhnite geometry if we let L —)■ cx) and set ipsioo) = 0. 

We next study the spectral symmetry as in the previous case. For the spin-hip operation, 
S, we hnd ^±[iSV’] = Thus, using the fact that u is unitary, the boundary 

condition is mapped to with u' = —A boundary condition preserving 

the spectral symmetry is found when u = u'. This latter condition can be rewritten as 
M = = M-^ for M = i(TzU. That is, M. is Hermitian and squares to 1. There 

are two types of solutions to this equation; (1) M. = ±1, yielding u = and (2) 

M. = rrixCTx + nTyCTj, -f rrizCTz = m • cr, with a unit vector m = {rrix, my, mz) G on the real 
two-sphere, yielding u = —iaz{m ■ cr). It is not difficult to show that in the hrst case, the 
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(43) 


boundary conditions for the wavefunctions read 

V’t(O) = ±'0^(0), = ±'ipi{L). 

As an example of the second case, we take u = =Fil. Then, 

V^t(O) = ±^|(0), = tA{L). (44) 

The effect on the spectrum can be seen directly in terms of the number of zero-energy 


states in each case. With the boundary conditions in Eq. (43) there is always one exact 
zero-energy state. 




(+), 


_ —mxlv 


X = e 


or 


- 1 . 


(45) 


For the the boundary conditions in Eq. (44) there are no exact zero-energy states for hnite L. 


However, for large L the conditions with the upper sign are satished by both and ipQ ^ 
to the order while for the lower sign the conditions are never satished independent 

of L. Thus, in this case, there are no zero-energy states for the lower sign and there are two 
asymptotically exact zero-energy states for the upper sign. 

Thus, we see that even after imposing the spectral symmetry given by S, the bulk¬ 
boundary correspondence may still be incomplete. This can be understood in terms of the 


physical implications of the boundary conditions in Eqs. (43) and (44). These conditions 


can be understood as modeling the interface of the wire with an environment. Note that the 
condition '?/’'|-(0) = ' 04 .( 0 ) fh® semi-inhnite case allowed the existence of a zero-energy state 
localized at the edge, while the condition 0'f(O) = —"0^(0) did not. Therefore, the former 
condition models the interface of a normal insulator, the “vacuum,” with the left edge of the 
wire that supports a topological bound state. If we had chosen the semi-inhnite geometry 
a; < 0 with a right edge, it is easy to see that the sign of the boundary condition modeling 
this situation would have been the opposite. This shows up in the hnite geometry: now the 


boundary conditions in Eqs. (43) and (44) at the two edges separately model an interface 


with vacuum with the -|- sign on the left edge and the — sign at the right edge. 


VII. SELF-ADJOINT EXTENSIONS AND CONSERVED CURRENTS 


As we have seen in several examples, such as in Eq. (42), the conditions determining the 


domain of self-adjoint extensions of a Hermitian operator were written in a way that was 
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independent of the specific choice of lA. In this section we show that this is quite a general 
result. Furthermore, we show that this can be directly interpreted as the conservation of a 
local current associated with the self-adjoint operator. 

Given a self-adjoint operator A and ^ we must have {'ip\A(l)) — = 0. 

This equation is the basis of the self-adjoint extensions we have so far discussed. For example, 
for the momentum operator P = —id/dx on the Hilbert space .^^[0,L] we have 


0 = (v^|P0) - {pm 


d 


—i—(j){x) — —i—'ilj{x) (j){x) 


= —i 


dx 

{4)* {x)(j){x\] dx. 
1 da; ^ 


d 


dx 


dx 


(46) 


As a result, we hnd '0*(T)0(L) —'0*(O)0(O) = 0. The general solution independent of -0 and 


0 is '0(T) = e*^'0(O) as we found under Eq. (25). Note that for -0 = 0 this condition reads 

|0(L)r = 10(0)10 (47) 


which is simply the conservation of probability at the two boundaries. Since the momentum 
operator generates space translations, the probability density can be interpreted as the 
conserved current under space translations. Therefore, this condition simply states that the 
probability density must be conserved throughout the interval [0,L]. 

Indeed, we show now that this way of specifying self-adjoint extensions of a Hermitian 
operator is quite general. For a self-adjoint operator A one can always dehne a one-parameter 
family of unitary operators Ua_{c() = generated by A. This result is known as Stone’s 

theorem .11^ For example, the Hamiltonian generates the time-evolution operator with a being 
the time parameter. Then we can “evolve” any given state 0 G ^(A) as 0(a) = Ua{oi)'iP 
with the initial condition 0(0) = 0. Since for two such evolved states, 0(a) and 0(a), the 
inner product (0(a) |0(a)) = / 0(a;, a)*0(a;, a) dx is independent of a, the density pA^x., a) = 
0(a;, a)*0(x, a) must have a local conserved current, jA{x,a), associated with it such that 
the following continuity equation is satished: 


d d 

jtPa + IpdA = 0 . 

oa ox 

To see this, note that since id'ilj{Q)/da = Atl^l^a), we have 
d 

i-A—pA{x, a) = ip*{x, a){A(j))(x, a) — (A0)*(a;, a)0(a;, a). 
oa 


(48) 


(49) 
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Thus, the self-adjoint condition can be directly written as 


0 = - {Aij\(j))] 

= - J -^pA{x,a) dx 

= f -^jAix,a)dx. 


( 50 ) 


Since this is satished for all a, and in particular a = 0, we hnd that the conserved current 


j^(x) =i [i’*{x'){A(j)){x') — (Ai/j)*{x')(j){x')] dx' 


(51) 


characterizes the domain of the self-adjoint operator A. 


For the continuum Hamiltonian "H, Eq. (12), one hnds the local current 


jn(x) = -tn/i(a:)V,<^(x) = iv[4!AA4>l(x) - i/>J(a:)<#'t(a:)] ■ 


(52) 


So, in the semi-inhnite geometry, noting that = 0, we hnd ^-^(0) = 0. This reads 


V’t(O)0t(O) = V’j)(O)0t(O)) which indeed yields Eq. (30) independent of the states. In the 


hnite geometry x G [0, L] we have jni.L) = ^-^(0), which is precisely the condition we found 


in Eq. (42). In order to obtain the boundary conditions from the current condition, we 


must revert the steps that took us from Eq. (41) to Eq. (42). To do so, we use the identity 


jn = + #t) - ('0t + + #t)] to write jn{L) = jn{0) as 

-V'^(O)+#t(0)y+#t(0)\ ^ V’i(O)+#t(0) y ^ 0i(O)+*0t(O) ^ 


(53) 


This means that the boundary conditions must map 


-V'i(o) + #t(o)^ ^ 

>^(^) + #t(^) / 


V’i(o) + #t(o) 

-^t(^) + #1(^)2 


(54) 


such that the inner product in Eq. (53) is preserved independently of the state 'ifj. This is 


satished if and only if the states satisfy the boundary condition in Eq. (41). This method 


can be used to obtain the boundary conditions from the conserved current more generally; 
see the discussion in Appendix [Cj 
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VIII. CONCLUSION 


As elaborated in this paper, the study of the continuum limit of a lattice model is a 
powerful method for extracting the universal, long-distance physics of the system. Different 
lattice models can share the same continuum, or held-theory, description and, therefore, the 
same universal behavior. In the context of topological phases, it is instructive to compare 
our results to those for the Kitaev model for a one-dimensional p-wave superconductor 
This is a model proposed for the realization of exotic topological bound states, known as Ma- 
jorana fermions, with special non-Abelian braiding properties potentially useful for quantum 
computation. It has recently received great attention as experimental efforts to realize the 
model in nanowires have produced encouraging results.!^ While the physics of this system 
is quite different from the one we studied in this paper (e.g., there is no superconductivity 
in our system), the continuum limits of the two systems are equivalent. This means they 
have the same topological phase diagram. In particular, the self-adjoint extensions in the 
continuum and the topological bound states are in one-to-one correspondence to each other. 
What differs from one continuum model to the other is the specihc structure of the quantum 
helds and, most importantly, the physical properties of the different phases. In fact, it turns 
out that these two models are not only equivalent in the continuum limit but that their 
lattice versions are unitarily equivalent (see Appendix for a proof). So, our lattice and 
held theory analyses are directly applicable to the study of topological phases of the Kitaev 
model. 

Topological bound states, such as Majorana fermions, are interesting in part for their 
usefulness in device applications. These applications are often based on dynamical features of 
the bound states as they are moved around or otherwise manipulated. In such manipulations, 
the Hamiltonian of the system is changed in time by varying the external parameters of the 
system. For example, by changing the Zeeman energy /i one can manipulate the position 
of the bound states. It is worth emphasizing that, when analyzing such schemes in the 
continuum, it is important to work with a self-adjoint extension of the Hamiltonian since 
only such an operator can dehne the time evolution of the system properly. 

In summary, through a simple model we have illustrated the importance of self-adjoint 
extensions in the continuum, or held-theory, description of topological phases of quantum 
systems with boundaries. In particular, we clarihed the physical interpretation of the ex- 
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tended operators in terms of a conserved local cnrrent. These extensions correspond to 
different physical situations with physically distinct environments outside the system, or 
equivalently, experimental conditions. We showed that the distinction can persist even after 
imposing internal symmetries that restrict the choice of the extension. Thus, the notion of 
bulk-boundary correspondence in a topological phase cannot in general be dehned indepen¬ 
dently from the choice of the extension and the corresponding boundary conditions at the 
edges of the system. 
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Appendix A: Continuum ladder operators 


Here we show that the operator 4/^ dehned in Eq. (14) is indeed a ladder operator 
satisfying the commutation relation [H^, Using the algebraic relation [AB, C] = 


A{BC -\r CB) — {AC -1- CA)B and the anticommutation relations (^, we have 


1^0. i-t 


dx dx' 


SlS2S' 

SlS2S' 

^ jH sis24>s'e{x')^1^{x)6s2sA{x - x')dxdx' 

SlS2S' 

/ 'Hsis'(I)s'e{x)^1^{x) dx 

sis' 

/ E(j)s,E{x)^l^{x) dx 


dx dx' 


Si 




(Al) 


where, in the penultimate line, we also used Eq. (13). 
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Appendix B: Dense and closed operators 


In this appendix we give definitions for densely defined and closed operators. We illustrate 
these concepts with the momentum operator P. 

Definition 1. A subset C is dense if for every ip G there is a 
sequence ipn & ^ that converges to ip in norm, written simply ipn —t ip, i.e., 
lim„^oo W'ipn - 'ipW = 0. 

Definition 2. An operator A : ^(A) —)■ is a densely defined operator if 

^{A) is dense in 

Equivalently, an operator is densely defined if there is an orthonormal basis (pj G 3>{A) 
for the Hilbert space such that any state ip G Jlf can be written uniquely as a superposition 
of (pj, that is. Ip = '^jCij(pj, with Oj = {pj\ip). Here, the sum is formally understood as 
a sum over discrete values of j and an integral over continuous values of j. 

An operator A must be densely defined in order to uniquely define its adjoint, Afi with 
the relation (plA'^ip) = {Ap\ip)] when this property holds for all pj, A^p is uniquely written 
as A^ = 

Definition 3. The operator A is closed if for any sequence pn G t^(A) for which 
Pn ^ P and Apn —t p, we have p G ^(A) and p = Ap. 

Let us show that the momentum operator P = —id/dx, with 

^{P) = {p & .if^[0,L] : p is absolutely continuous,■^(0) = p{L) = 0}, (Bl) 

is densely defined and closed. First note that there is an orthonormal basis pj G 3’{P), 
j G N, given by pj{x) = a/ 2/L sin(j7ra;/L). For any p G Jif, the sequence of partial Fourier 
series pn defined by 

n pL 

Pn{x) = 'S^ajpj{x) with Oj = / p*{x)p{x)dx (B2) 

i=i 

converges to p in norm. Thus, ^{P) is dense in .if^[0,L]. 

To show that P is closed, we assume a sequence Xn G ^{P) is given so that Xn ^ P 
and Pxn —t p- Since ^{P) is dense, we may write p = Ej and define the partial-sum 


22 



sequence V’n = YTj ^ ^(-P) that also converges to tjj, V’n V'- By continuity, Pijjn —)■ 0 , 
too. Since lim„_,.oo P’^n = exists, lim,^^oo 0n = 0 is absolutely 

continuous and vanishes at the boundaries, a: = 0, lP therefore, V’n —)■ 0 G ^{P). Finally, 
we show that P0„ —)■ Pip and, thus, 0 = Pip: 



i’n{x) 




Jtvx 

~T~ 



j sm ——ip[x ) ax 


cos 


JTTX 


d jnx' 
— cos — 


ip{x') dx' 


cos 


JTIX 


J7TX JTTX 

COS —— cos —— 
1j Jj 


d 


— 0 (a;') dx'- 


since 



i=i 


we conclude that 

lim ||P0n — P0II = 0. 

n^oo 


(B3) 


(B4) 


(B5) 


Appendix C: Boundary conditions from conserved current condition 


In order to find the boundary conditions from the conserved current of Sec. VII , we note 
that the current Ja is a local, sesquilinear, Hermitian form of the two states 0, 0, such that 
jAi'ipyfp) = 7 ^( 0 , 0 )*- A general expression for such a form is Ja^x) = i/[ 0 (a;)]l Ja i/[ 0 (x)], 
where u[ip{x)] is an n-dimensional vector of linear combinations of the elements of ip{x) 
and its derivatives, and j\ = J 4 is a Hermitian matrix. For example, for the momentum 
operator P = —id/dx, we have i/[0(a;)] = ip{x), n = 1, Jp = 1. 


Exercise 2. For the operator A = —dJ/dx'^, show that z/[0] = (0 dip/dxY, 

Ja = CTy, and n = 2. 

We can diagonalize Ja = T^PT, with T unitary and D diagonal with real elements. The 
unitary matrix T can be always chosen such that the diagonal elements of D are sorted in 
decreasing order. Then D = (d+ © d_)^, with © the direct sum, such that d± = are 
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Hermitian and anti-Hermitian, respectively, and diagonal; the diagonal elements of d± are 
simply the square roots of, respectively, the positive and negative elements of D. 


Exercise 3. For the continuum Hamiltonian, Ti = ivayd/dx + show that 
h'l'ip] = 'ip, = —vay, and n = 2. Also, show that T = = (1 + 


D = vaz = 


V 0 


(Cl) 


0 —iy/vj 

and, thus, dy = a/F, d_ = —iy/v (up to signs). 

Now, writing Ja = T'^idy © d_)‘^T = [{dy © (i_)'*'T] ^ [((i+ © d-)T], we can see that 

Ja = {{d+ © -d-)T i/[^/>])^((d+ © d-)T i/[0]) 

= (C2) 

where we have dehned 

{d+® d_)T ulpj] =: Uylpj] ® (C3) 

We will now consider two different geometries. First, consider a semi-inhnite geometry 
X G [0, cx)). Then the current condition jA(0) = 0 reads 

z/+[?/>(O)]'^z/+[0(O)] = z/_[?/>(O)]V_[0(O)]. (C4) 


If, and only if, dy and d- (and, thus, z/+ and i/_) have the same dimension, n/2, this 
condition has a solution in terms of a unitary f x | matrix u that maps h->■ z/_, i.e.. 


uiy+['i/j{0)] = z/_[V'(0)]. 


(C5) 


Otherwise, A will have no self-adjoint extensions. For example, for the momentum operator 
P = —id/dx, we have z/+['F’(0)] = V’(O) s-nd z/_ does not exist. Therefore, P cannot have a 
self-adjoint extension on the semi-inhnite line. For the continuum Hamiltonian, Ti, we have 
z/±['0(O)] = y/vj2 ['0t(O) i *'04.(0)]- Therefore, the self-adjoint extensions of "H are character¬ 


ized by a U(l) phase u = after simple algebra, Eq. (C5) reads '0-f(O) = cot(6'/2)'F’4,(0), 


which is the same as Eq. (30) with 6 = 7r/2. 

Now consider a hnite geometry x G [0, Lj. The current condition, Ja^L) = 4 ( 4 ( 0 ), yields 

z/+[?F(O)]V+[0(O)] + z/_[V'(L)]^z/_[(F(L)] = z/_[?F(0)]^z/_[(F(0)] + z/+[?F(L)]'^i/+[0(L)] (C6) 




(C7) 
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where 


= i"±[V'(0)] © z/^['0(L)]. (C8) 

Thus, the state-independent boundary conditions are given by a norm-preserving, n x n 
unitary matrix u that maps h- i.e.. 




(C9) 


In the case of the momentum operator P = —id/dx we see immediately that now <^+["0] = 
-0(0) and Thus, the self-adjoint extensions of P are characterized by a U(l) 

phase as iIj{L) = e*^'0(O)- the case of the continuum formulation of the topological wire, 
we hnd 


‘^+ = -i 






^t(0)+#t(0) ^ 


Thus, from Eq. (C9) we hnd the boundary conditions in Eq. (41). 


(CIO) 


Appendix D: Unitary equivalence between a ferromagnetic insulator and a p-wave 
superconductor 


Here we prove that the ferromagnetic insulator model in Eq. ([^ is unitarily equivalent 
to a p-wave superconductor. More precisely, the superconducting model represents two 
uncoupled Kitaev chaind^ that share identical energy spectra. 

Dehne the local canonical map 


t _ /r + idl \ _ fr + Wr 

V2 ' V2 ' 


(Dl) 


in terms of two new sets of fermion creation and annihilation operators, /I and /^, and gl and 
Qj.. It is straightforward to show that the / operators commute with the g operators and each 
pair separately satishes fermionic anti-commutation relations {/^, /J,} = {g^, gl,} = drr', etc. 
Using these relation and after some algebraic manipulations, the Hamiltonian H of Eq. ([^ 
becomes 


N r 


h + ,,n^Y1 


r=l 


^{Prfr + aldr) + Hflfr+l + aldr+l) + fl+1 + dlgl+l) 


X 


-hh.c., (D2) 
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where, as before, the sum is understood not to contain any terms with the index + 1. 
Thus, the Hamiltonian is mapped to two copies of the Kitaev model for / and g separately. 
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